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Abstract 

The  ancient  Babylonians  had  an  iterative  technique  for  numerically 
approximating  the  values  of  square  roots.  Their  method  can  be  physically 
implemented  using  series  and  parallel  resistor  networks.  A  recursive  formula 
for  the  equivalent  resistance  RQq  is  developed  and  converted  into  a  nonrecursive 
solution  for  circuits  using  geometrically  increasing  numbers  of  identical 
resistors.  As  an  example,  24  resistors  R  are  assembled  into  a  second-order 
network  and  RQq/R  is  measured  to  equal  y/l  to  better  than  0.2%,  as  could  be 
done  in  an  introductory  physics  laboratory. 

(Some  figures  may  appear  in  colour  only  in  the  online  journal) 


1-  Aim 


Given  a  bucket  of  identical  resistors  R ,  construct  a  circuit  whose  equivalent  resistance 
approximates  \flR  (to  any  desired  degree  of  accuracy,  limited  only  by  the  precision  of  the 
resistors).  To  ensure  that  the  construction  is  accessible  to  introductory  physics  students,  the 
resistors  are  to  be  combined  either  in  series  or  in  parallel  at  each  stage  of  the  circuit  [1]. 

2.  Theoretical  analysis 

To  simplify  the  equations  that  follow,  choose  units  of  resistance  so  that  R  =  1 .  In  the  next 
section,  it  will  be  explained  how  to  restore  conventional  units  of  ohms  to  the  terms. 

Start  by  putting  two  resistors  in  series  to  get  a  resistance  of  2.  Denote  that  combination  as 
So  where  the  ‘S’  refers  to  putting  two  initial  resistors  in  series  and  the  subscript  ‘0’  denotes  a 
zeroth-order  approximation  to  \fl.  However,  Sq  is  larger  than  \fl  and  is  thus  an  overestimate. 
On  the  other  hand,  one  gets  an  underestimate  by  constructing 


(1) 


where  in  later  steps,  ‘P’  will  refer  to  putting  two  previous  resistor  combinations  in  parallel. 
These  two  zeroth-order  circuits  are  illustrated  in  figure  1 . 
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50  =  2 


P0=l 

Figure  1 .  Starting  circuits  to  construct  a  resistor  network  whose  equivalent  resistance  approximates 

V2. 


^  =  3/2 


Pi  =  4/3 


Figure  2.  The  two  circuits  that  can  be  constructed  using  six  unit  resistors  whose  equivalent 
resistances  are  nearest  in  value  to  >/2. 


To  achieve  an  equivalent  resistance  that  is  nearer  in  value  to  \/2,  one  can  average  together 
So  and  Po  to  get 

Si  =  ^(S0+P0)  =  l-  (2) 

To  accomplish  that,  first  put  So  and  Po  in  series  with  each  other.  Then  put  that  trio  of  unit 
resistors  in  parallel  with  an  identical  trio  to  halve  the  total  resistance.  On  the  other  hand,  if 
one  first  puts  So  and  Po  in  parallel,  and  then  puts  that  trio  in  series  with  an  identical  trio,  one 
doubles  the  resistance  of  each  trio  to  end  up  with 


p1=2fi+iy1=2^- 

Vs0  Po)  s0  +  p0 


SqPq 

~sT 


(3) 


using  equation  (2)  in  the  last  step.  Consequently 


S\P  i  =  S0Po  =  2 


(4) 


from  equation  (1)  and  therefore 


(5) 


These  two  first-order  resistor  combinations  are  sketched  in  figure  2.  Again  S i  is  an 
overestimate  and  Pi  is  an  underestimate  to  the  value  of  yfl. 
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Table  1 .  Values  predicted  by  equation  (9)  for  a  starting  value  of  So  =  2. 


n 

Sn 

0 

2 

1 

3/2=  1.5 

2 

17/12  «  1.417 

3 

577/408  «  1.414  22 

These  two  constructions  can  now  be  iterated  to  successively  approach  an  equivalent 
resistance  of  V2.  To  form  S2  and  P2 ,  simply  replace  every  occurrence  of  So  and  Po  in  figure  2 
with  Si  and  P\,  respectively.  Then  equation  (2)  becomes 

fc  =  i(S.+fl,  =  i(s.  +  £)  (6) 

using  equation  (5)  in  the  last  step,  and  equation  (4)  becomes 

S2P2  =  S1P1  =  2  (7) 

so  that 

P2  =  2  (8) 

S2 

It  takes  24  unit  resistors  to  build  either  of  the  second-order  networks  S2  or  P2. 

More  generally  for  n  ^  1  one  can  use  1.5  x  4n  unit  resistors  to  construct  the  two  circuits 


and 


(9) 


Pn=^-  (10) 

Equivalent  resistances  computed  recursively  using  equation  (9)  are  listed  in  table  1 .  One 
sees  that  Sn  rapidly  converges  to  ~  1.41421. 

Recursion  relation  (9)  can  be  converted  into  a  nonrecursive  formula  valid  for  n  ^  0, 


s„  =  V2 


l  _|_  rf(n) 
l  _  rf(n ) 


(11) 


where  f(n)  =  2n  and  r  =  (V2  —  l)/(\/2  +  1)  ~  0.1716,  as  verified  in  appendix  A.  It  follows 
immediately  from  this  formula  that  Sn  —>  \fl  as  n  00.  One  can  more  compactly  write  this 
result  as 


Sn  =  V2  coth  (2"acothV2)  (12) 

where  coth  and  acoth  are  the  forward  and  inverse  hyperbolic  cotangent  functions.  (This  formula 
again  implies  Sn  —>  \fl  as  n  oo  because  coth  oo  =  1 .)  This  explicit  formula  for  Sn  is  plotted 
in  figure  3  as  the  solid  curve,  treating  n  as  a  continuous  variable.  Its  values  for  the  integers 
n  =  0  to  5  are  indicated  by  the  dots.  The  asymptotic  value  \fl  is  denoted  by  the  horizontal 
dashed  line. 

To  prove  that  equation  (9)  converges  to  without  using  equations  (1 1)  or  (12),  assume 
[2]  that  a  limiting  value  of  Sn-  \  exists  and  call  it  S.  That  must  also  be  the  limiting  value  of  Sn 
and  hence  equation  (9)  in  the  limit  n  — >►  oo  becomes 


S 


1 

2 


=>•  S2  =  2 


(13) 
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n 

Figure  3.  Explicit  and  recursive  values  of  Sn  asymptotically  approaching  the  value  a/2. 


whose  positive  solution  is  S  =  (For  another  approach,  see  appendix  B.)  This  limiting 
value  can  be  understood  intuitively  by  writing  equation  (9)  as 

Sn=l-{Sn^+Pn-X).  (14) 

Since  Sn-\  is  slightly  larger  than  a/2,  it  follows  that  Pn_\  =  2/ Sn-\  will  be  slightly  smaller  than 
a/2.  Therefore  averaging  together  these  two  values  gives  an  improved  estimate  for  a/2.  This 
iterative  method  of  approximating  a/2  was  known  to  the  ancient  Babylonians  [3],  motivating 
the  title  of  this  paper. 

3.  Experimental  verification 

Students  can  assemble  one  of  these  Babylonian  resistor  circuits  and  measure  its  resistance. 
Noting  from  figure  3  that  good  accuracy  is  already  obtained  for  n  =  2  and  that  it  would  be 
difficult  to  correctly  wire  96  or  more  resistors  together,  it  is  reasonable  to  build  a  second-order 
network.  Since  Pn  is  closer  in  value  to  a/2  than  is  Sn  (for  any  value  of  n ),  it  makes  sense 
to  construct  P 2  rather  than  S2.  All  of  the  equations  in  section  2  assumed  unit  resistors;  for 
laboratory  resistors  of  resistance  R ,  one  simply  has  to  multiply  any  value  or  expression  for  Sn 
or  Pn  by  R.  (For  example,  S 1  in  figure  2  then  has  a  resistance  of  1.5 R.) 

The  required  components  are  24  identical  resistors  R  and  five  short  jumper  wires.  First 
wire  the  resistors  together  as  four  parallel  lines  of  6 R  each  by  twisting  their  ends  together 
(and  securing  them  with  solder),  and  then  add  the  jumper  wires  indicated  in  red  in  figure  4. 
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Figure  4.  Circuit  P 2  redrawn  to  make  it  easy  to  assemble. 


Figure  5.  Measurement  of  the  equivalent  resistance  of  circuit  P 2. 


The  equivalent  resistance  between  the  left  and  right  ends  is  predicted  to  be 
24  r~ 

—R  «  V2R.  (15) 

The  difference  between  24/17  and  y/2  is  less  than  0.2%. 

The  circuit  in  figure  4  was  constructed  using  24  nominally  1  kQ  resistors.  The  actual 
resistances  were  measured  using  a  multimeter  prior  to  assembly  and  were  found  to  be 
R  =  983  zb  5Q.  Thus,  either  the  resistors  are  systematically  low  or  the  multimeter  is  a 
bit  out  of  calibration,  but  that  will  not  affect  the  ratio  RQq/R  as  long  as  the  same  multimeter  is 
used  to  measure  the  equivalent  resistance  RQq.  As  can  be  seen  in  the  photograph  of  the  setup 
in  figure  5,  we  measured  Req  =  1 388^2,  which  differs  by  less  than  0.2%  from  \fl R  =  1390  £2 
and  is  well  within  the  5  Q  variation  in  the  individual  resistors. 

4.  Conclusions 

Discussion  and  assembly  of  these  resistor  networks  provide  an  interesting  alternative  to  routine 
series  and  parallel  circuits  in  an  introductory  physics  course.  An  extension  of  the  results 
presented  here  is  to  keep  P0  =  R  (in  laboratory  units)  but  replace  So  with  any  resistance  kR 
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instead  of  2 R.  The  equivalent  resistance  of  networks  iterated  in  the  style  of  figure  2  will  then 
asymptotically  approach  the  value  VkR  rather  than  \/2 R.  The  number  k  need  not  even  be  an 
integer  but  can  have  any  real,  positive  value.  (For  example,  if  Po  =  200£2  and  So  =  4.7  k£2, 
then  k  =  Sq/Pq  =  23.5.)  Such  Babylonian  networks  can  thus  be  thought  of  as  an  analogue 
calculator  of  square  roots.  Another  option  is  to  replace  the  resistors  with  capacitors  or  springs 
because  they  too  obey  series  and  parallel  combination  rules. 


Appendix  A.  Verification  of  equations  (11)  and  (12) 


Since  f(n  +  1)  =  2  •  2n,  equation  (11)  implies  that 


Sn+ 1  =  >/2 


1  +  r2f 

1^7' 


(A.l) 


But  according  to  equations  (9)  and  (11)  this  result  is  supposed  to  be  equal  to 

1  /  2  \  _  t  1  +  Y  1  1  -  Y 

2  V” +  s~n)  -  7!  +  vf  T+7r 

In  other  words,  one  needs  to  show  that 


(A.2) 


^  1  +r2f  _  1  (1  +  rf)2  1  (1  -  rf)2 
V2  1  -r2/  “  71  1  -r2/  +  71  1  -  r2f  ' 


(A. 3) 


But  this  equality  is  verified  by  expanding  the  two  squares  in  the  numerators  of  the  right- 
hand  side.  Different  values  of  r  give  rise  to  different  initial  values  for  Sq.  In  particular, 
r  =  (\/2  —  l)/(\/2  +  1)  results  in  So  =  2. 

To  directly  obtain  equation  (12),  run  the  recursion  solver  command 


RSolve[{S[n  +  1]  =  =  S[n]/2  +  1/S[n],  S[0]  ==  2},  S[n],  n]  (A.3) 


in  Mathematica.  Alternatively  one  can  convert  between  equations  (11)  and  (12)  using  standard 
exponential  and  logarithmic  identities  for  coth  and  acoth. 


Appendix  B.  A  physics-based  approach  to  finding  the  limiting  value  of  Sn 


Rewrite  equation  (9)  as 

‘Sft+l  _  1  Sn 

(w+  1)  -n  ~  ~Sn  ~  ~2m  {  ' 

The  left-hand  side  is  ASn/An  ~  dSn/dn  by  treating  Sn  and  n  as  continuous  variables. 
Substituting  that  derivative  into  equation  (B.l)  and  multiplying  through  by  2 Sn,  one  obtains 


dOQ 

dn 


2 -St 


(B.2) 


By  viewing  Sn  as  a  velocity,  equation  (B.2)  can  be  interpreted  as  the  change  in  kinetic  energy 
of  a  particle  falling  in  a  uniform  gravitational  field  while  subject  to  a  quadratic  drag  force.  The 
differential  equation  is  separated  and  integrated  to  get 


Sn  as  a/2  -  e-". 


(B.3) 


The  limit  n  — >■  oo  corresponds  to  a  terminal  speed  of  a/2  in  which  case  the  left-hand  side  of 
equation  (B.2)  is  zero  since  the  velocity  is  no  longer  changing.  Equation  (B.3),  unlike  equations 
(11)  and  (12),  is  only  approximate  for  finite  n ,  however.  In  contrast  to  the  red  curve  in  figure  3, 
equation  (B.3)  describes  a  graph  which  starts  at  1  and  increases  towards  an  asymptotic  value 
of  >/2. 
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